Fermi surface and heavy masses for UPd2Al3 
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We calculate the Fermi surface and the anisotropic heavy masses of UPd2Al3 by keeping two 
of the 5/ electrons as localized. Good agreement with experiments is found. The theory contains 
essentially no adjustable parameter except for a small shift of the position of the Fermi energy of 
the order of a few meV. A discussion is given why localization of two / electrons is justified. 



There has been experimental evidence since a number 
of years that in UPd2Al3 the 5/ electrons have a dual 
character. The large observed magnetic moment of 0.85 
/is [1] suggests localization of the / electrons. On the 
other hand, the large jump in the specific heat at the 
superconducting transition temperature Tc=2K implies 
a large Sommerfeld specific heat coefficient 7 and mod- 
erately heavy fermion behavior [2]. This requires delo- 
calized 5/ electrons. The dual character is also found 
in a number of other physical properties such as pho- 
toemission, inelastic neutron scattering and muon spin 
rotation [3-6]. It is also observed in other metallic U 
compounds [7,8]. Moreover, the assumption is supported 
by quantum-chemical calculations on U(C8H8)2 [9] which 
exhibit a number of low- lying excitations caused by intra- 
atomic rearrangements of the 5/ electrons. There are 
speculations that similar 5/ excitations might even be 
responsible for the formation of Cooper pairs in supercon- 
ductors [10]. In that case the dual model should allow for 
a natural description of heavy-fermion superconductivity 
coexisting with magnetism produced by 5/ electrons. 

We have previously applied the dual model in order 
to explain the mass enhancement of the quasiparticle 
excitations and the de Haas-van Alphen frequencies in 
UPta, the model system of heavy-fermion behavior in U 
compounds [11]. Our theory conjectured that the de- 
localized 5/ states hybridize with the conduction elec- 
tron states and form energy bands while the localized 
ones form multiplets in order to reduce the local mutual 
Coulomb repulsions. The interaction of the two subsys- 
tems, i.e., the delocalizcd with the localized one leads to 
a mass enhancement of the former. The situation resem- 
bles that of Pr metal where a mass enhancement of the 
conduction electrons by a factor of 5 results from virtual 
crystal- field excitations of the localized 4/ electrons [12]. 
The same mechanism leads to the heavy quasiparticles 
in the recently discovered heavy-fermion superconductor 



PrOs4Sbi2 [13]. The role of the incompletely filled 4/ 
shell is taken by an incompletely filled subshell of 5/ 
electrons when U ions are considered instead. 

The coexistence of itinerant and localized 5/ states 
is referred to as partial localization. It plays an impor- 
tant role in many intermetallic actinide compounds. Par- 
tial localization arises from an interplay between the hy- 
bridization of the 5/ states with the environment and 
on-site Coulomb correlations. This is discussed below 
in more detail. The underlying microscopic origin is an 
area of active current research [14-16]. The associated 
physical picture is therefore quite different from the one 
suggested in Ref. [2], where two disjunct subsystems in 
k-space were postulated. 

The dual model should be contrasted with density 
functional based calculations which use the local den- 
sity approximation. These have been successful in ex- 
plaining the measured de Haas-van Alphen frequencies 
of systems like UPd2Al3 or UPta [17-19], but they fail 
to explain the observed heavy quasiparticle masses. For 
UPta the observed masses are by a factor of 20 larger 
than the calculated ones and for UPd2Al3 the difference 
is roughly a factor of 4. When using the self-interaction 
corrected local spin-density approximation (SIC-LSDA) 
[20] a ground state is found with coexisting localized (f^) 
and delocalized 5/ electrons of U. But the calculated den- 
sity of states is too small by a factor of approximately 10 
to account for the observed linear specific heat. 

While the calculated energy bands are too broad for ex- 
plaining the effective masses, they are too small in order 
to fit the observed photoemission data of UPta, UPd2Ala 
or UBcia [21]. The latter shows a broad peak just be- 
low the Fermi energy Ep and is quite different from the 
data of heavy quasiparticle systems which involve Ce'^+ 
instead of U ions. Here a broad peak is found approx- 
imately 2eV below Ep [22] while at the Fermi energy a 
small Kondo resonance is detected. We add that even for 
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UPda which has locahzed 5/ electrons and does not show 
heavy-fermion behavior, the photoemission data resem- 
ble that of the previously mentioned U compounds, ex- 
cept that there is no / weight at Ep. This indicates that 
despite the afore mentioned successes the strong correla- 
tions in those materials are not properly treated by the 
presently used computational methods. Instead we find 
that a more microscopic understanding of heavy-fermion 
behavior is highly desirable. 

The aim of the present investigation is to show that 
the de Haas-van Alphen frequencies for the heavy- 
quasiparticle portion of the Fermi surface and the large 
effective masses (including their anisotropics), can be ex- 
plained very well by treating two of the 5/ electrons as 
being localized. We put the localized electrons into 5/ 



j=| orbitals with j 
are treated as itinerant electrons 



±1 and ±i. 



Thej, 



±^ states 



The reason for this 
choice is explained below, but it is worth pointing out 
here that the jz = ±| states hybridize strongest among 
the different ones in a conventional LDA calculation. 

The calculations of the heavy bands proceed in three 
steps. First, the bandstructure is determined by solving 
the Dirac equation for the self-consistent LDA potentials, 
thereby excluding the U 5/ states with j=|, iz — ±| 



and jz = ±2 from forming bands. Two 5/ electrons in 
localized orbitals are accounted for in the self-consistent 
density and, concomitantly, in the potential seen by the 
conduction electrons. The intrinsic bandwidth of the itin- 
erant U 5/ j=f , jz = ±1 electrons is taken from the LDA 
calculation while the position of the corresponding band 
center C is chosen such that the density distribution of 
the conduction electrons remains the same as within the 
LDA. The / electron count per U atom for the delocal- 
ized 5/ electrons amount to Uf ~ |, i.e., the system is of a 
mixed valence type. The calculated de Haas-van Alphen 
frequencies are shown in Fig. 1 which also contains the 
experimental results. 
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FIG. 1. Comparison of the measured [19] (black symbols) 
and calculated de Haas- van Alphen frequencies for UPd2Al3. 
The dominant part reflects a nearly cylindrical sheet of the 
Fermi surface. 

One notices that the agreement for the heavy quasipar- 
ticle branches is very good. The frequencies referring to 
the light parts of the Fermi surface are less well repro- 
duced, but that is of no surprise. One cannot expect that 
die LDA reproduces well the relative shifts of the centers 
of the light bands which strongly effect the shape of the 
corresponding parts of the Fermi surface. The Fermi sur- 
face of the heavy quasiparticles in the antiferromagnetic 
phase consists of three sheets two of which are displayed 
in Figure 2. 
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We assign the orbit a to the first sheet (K-centered el- 
hpsoid) not shown in Fig. 2. The third sheet is the most 
important one and consists of a corrugated cyhnder and 
a H-centered ellipsoid (see Fig. 2a) to which we assign 
the orbits 7 and (3 and 62, €3, respectively. The branch 
C, is assigned to the second sheet whose shape depends 
sensitively on the position of the Fermi energy (see Fig. 
2b). For example a shift of the latter by 40 K changes 
that sheet from the one shown in Fig. 2b to that of Fig. 
2c. But the corresponding changes in the effective mass 
remain small. 

In the second step we calculate the multiplet struc- 
ture of the localized p states. This is done within 
the jj-coupling scheme, because the spin-orbit splitting 
is rather large. Thus a 6 x 6 Coulomb matrix has to 
be diagonalized for the two-particle states built from 
\j = |, U = ±1) and |j = |, = ±i) for the 
/2 subshell. The Coulomb matrix elements are calcu- 
lated following Condon and Shortley. Inputs are the 
Slater-Condon parameters F^ (Coulomb integrals) and 
(exchange integrals). The latter are evaluated with 

the radial function 5 (r) for U as obtained from a self- 
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consistent bandstructure potential. Thereby the chosen 
energy is that of the center of gravity of the 5/ bands. 
The required integrations are done within the atomic 
sphere surrounding the U ion. We note that we use the 
same Coulomb matrix as for UPta [11] where also two 5 
/ electrons with = ±|, ± ^ are considered. This is 
reasonable in view of the fact that in the nonrelativistic 
case the matrix elements agree up to two decimals with 
the ones computed in [9] for U(C8H8)2, indicating that 
they are insensitive to the chemical environment of U. 
The resulting eigenstates of the Coulomb matrix are no 
longer eigenstates of the total angular momentum J, but 
remain eigenstates of J^. 

We find a doubly degenerate ground state with J2=±3. 
It must be an cigenstate to J=4 since the Paiili principle 
requires an even value of J and J=0,2 are excluded. The 
states li = f , i = |, = 4, Jz = ±3) have an over- 
lap of 0.865 with the state ^H4. The latter follows from 
Hund's rule when the LS coupling scheme is applied. The 
two-fold degeneracy of the ground state is lifted by the 
crystalline electric field (CEF) which is acting on the 5/ 
subshell with the localized electrons. From experiments 
a r4 ground state has been suggested [23] so that the 
ground- and first excited state are 

I r3,4) = (I J = 4; J, = 3)± I J = 4; J, = -3)) . (1) 

For the splitting energy S between the two states a 
value of 20meV has been previously suggested [1] while 
a more recent value is 7meV [10]. The latter value will 
be used in the following. We want to point out that the 
next higher doublet of the Coulomb matrix is one with 
Jz = ±2 with an excitation energy of approximately 0.4 



cV. Therefore we may neglect all higher excited states 
and consider the two singlets ]r4) and [Fa) only. 

In a third and final step we determine the effective 
masses which result from the interaction of the delocal- 
ized 5/ electrons with the localized one. The renormal- 
ization of the band mass mb is given by 



m 



(2) 



w=0 



where denotes the local selfenergy of the delocalized 

5/ states. The latter is obtained by analytic continuation 
from the Matsubara frequencies = 7rT(2n-|-l) at the 
temperature T where it is given by 



in terms of the local susceptibility 



X{ien - iCn') 



-tanh 
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2T (^ie^-ie^,)2_s2 



and the local propagator 



I 



dE- 



N{E) 



le. 



-E- E{ien) 



(3) 



(4) 



(5) 



Here 2N(E) is the total density of states at the energy E 
as obtained from the bandstructure, when two 5/ elec- 
trons are kept localized. The prefactor a denotes the 5/ 
weight per spin and U atom of the conduction electron 
states near Ep. The matrix element M describes the tran- 
sition between the localized states )r4) and [Fa) due to 
the Coulomb interaction Ucoui with the delocalized 5/ 
electrons. It is given by 

M = if; 1^1® (r4 [C/couil Ts) ® 1/^ I ^) (6) 

and is directly obtained from the cxpection values of the 
Coulomb interaction in the 5/"^ states. The difference 

if; I, 11 ® iP; 4, 3 jL/couil 4, 3) ® \f\ |, |) - 
if; I, |l®(/'; 4, -3 jt/couil .P; 4, -3)®|/i; |, I) 
is -0.38eV. From this we obtain M=-0.19eV. When we 
want to write M in form of an exchange coupling be- 
tween the delocalized and localized 5/ states we need to 
know the Landc factor g. We find gcir=0.63 and with 
this value an exchange integral of size I~leV. This is of 
the correct size for 5/ electrons. 

In order to evaluate Eq. (2) we need to know N(E) 
and a. We extract the quantities from LDA calculations 
with two 5/ electrons being localized. For simplicity, we 
model N(E) by a Lorentzian. The actual value at the 
Fermi level N(0)=2.76 states /(eV cell spin) corresponds 
to the one found by [20]. A value of a2=0.44 is obtained 
and the value (5=7meV is used. 

The selfconsistent calculation yields a mass enhance- 
ment of 9.6. The resulting calculated quasiparticle 
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masses are in excellent agreement with experiment (see 
Table I). 



TABLE I. Effective masses for H 11 c. The experi- 
mental data are taken from Inada et al. [20] . 



m* (exp) m* (theory) 



C 65 59.6 

7 33 31.9 

13 19 25.1 

£2 18 17.4 

£3 12 13.4 

a 5.7 9.6 



What remains to be discussed is the justification for 
treating two of the 5/ electrons in orbitals jz = | and 
^ as localized. As pointed our before the hybridization 
of the jz = ±1 orbital with the neighboring atomic or- 
bitals is larger than the one for the orbitals with jz = ±| 
and So why are we allowed the neglect the hy- 

bridization of the latter orbitals altogether? The answer 
is found when the effects of intraatomic correlations are 
taken into account, i.e., those resulting from the on-site 
Coulomb and exchange effects. For a demonstration a 
two-site model was treated in Ref. [14] with anisotropic 
hybridization between the two U sites. An intermediate 
valency of 2.5 was used for the two atoms. Let \tpo) de- 
note the ground state of that system and ta the hopping 
parameters for the different orbitals. It was found in Ref. 
[14] that whenever one hopping parameter ta dominates 
the others, i.e., tee ^ tce',ta", the corresponding ground 
state expectation values is 

_ (^0 |t„.c+(l)c^.(2) IV^o) tg, 

(Vo|E^aC+(l)c„(2)|Vo) Eia ' 
a a 

i.e., the effective hybridization anisotropics are 
strongly exhanced by intra-atomic interaction. The op- 
erators c^{i){ca{i)) create (destroy) an / electron on site 
i in orbital jz = ol. The smaller hopping matrix elements 
are therefore suppressed. 
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